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STATISTICAL  FINE  STRUCTURE  IN  THE 
INHOMOGENEOUSLY  BROADENED  ELECTRONIC  ORIGIN 
OF  PENTACENE  IN  />-TERPHENYL 

*  .+ 

T.  P.  Carter  ,  M.  Manavi  ,  and  W.  E.  Moerner 
IBM  Almaden  Research  Center 
San  Jose,  California  95120 

ABSTRACT:  Recently,  the  first  observation  of  statistical  fine  structure  on  an 
inhomogeneously  broadened  absorption  profile  was  reported  [  W.  E.  Moerner  and  T.  P. 
Carter,  Phys.  Rev.  Lett.  59,  2705  (1987)  ]  for  mixed  crystals  of  pentacene  in  p-tcrphenyl 
using  laser  frequency-modulation  spectroscopy.  Statistical  fine  structure  is  time-independent 
structure  on  the  inhomogeneous  line  caused  by  statistical  variations  in  the  spectral  density 
of  absorbers  in  each  frequency  interval.  In  this  work,  a  model  and  an  analysis  of  statistical 
fine  structure  using  autocorrelation  techniques  are  presented,  and  the  dependence  of  the  effect 
for  pentacene  in  />-terphcnyl  at  1.4  K  on  modulating  frequency,  detection  phase,  center 
concentration,  and  position  in  the  sample  is  described.  Statistical  fine  structure  provides  a 
new  method  for  probing  inhomogeneous  systems  that  allows  investigation  of  the  underlying 
probability  distribution  function  and  the  determination  of  the  homogeneous  linewidth.  This 
fine  structure  also  represents  a  fundamental  limit  on  the  detectability  of  shallow  spectral 
features  in  inhomogeneously  broadened  lines.  ...  ....  — — 
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I.  Introduction 


Inhomogeneous  broadening  (IB)  of  atomic  or  molecular  transitions  occurs  when  the 
individual  absorbers  in  a  collection  possess  homogeneous  linewidths  which  are  narrower  than 
the  overall  lineshape  of  the  collection  as  a  whole,  due  to  a  distribution  of  individual  center 

frequencies.  IB  is  a  common  feature  of  the  low-temperature  spectra  of  impurity  centers  in 

2 

solids  that  appears  in  optical  spectroscopy  as  well  as  in  nmr,  esr,  and  Mossbauer 
experiments.  3  For  zero-phonon  electronic  transitions  of  crystal  impurities,  the 
inhomogeneous  distribution  (for  a  particular  impurity  orientation  or  site)  is  a  result  of  any 
of  a  variety  of  imperfections  in  the  host  crystal  such  as  dislocations,  point  defects,  or  random 
internal  electric  fields  and  field  gradients.  In  amorphous  hosts,  it  is  principally  the  large 
multiplicity  of  available  sites  that  leads  to  IB. 


IB  is  often  regarded  as  an  obstacle  to  the  determination  of  the  underlying  homogeneous 
width,  which  is  a  fundamental  property  of  the  impurity-host  combination  at  a  given 

4 

temperature.  Several  powerful  spectroscopic  methods  have  been  developed  to  extract  this 
information  from  inhomogeneous  systems,  including  persistent  5  and  transient  hole-burning 
as  well  as  such  coherent  transient  techniques  6  as  photon-echo,  optical  frce-induction  decay 
and  optical  nutation.  Fach  of  these  methods  has  its  strengths  and  weaknesses;  for  example, 
not  all  systems  exhibit  persistent  hole-burning,  and  some  that  do  are  complicated  by  spectral 
diffusion  mechanisms  which  prevent  the  direct  determination  of  the  homogeneous  width.  If 
hole-burning  docs  occur,  the  coherent  transient  techniques  may  not  be  applicable  because 
the  high  peak  intensities  normally  required  in  these  experiments  lead  to  rapid  hole-burning 
and  loss  of  signal.  However,  IB  is  not  always  regarded  in  a  negative  sense:  its  existence  is 
fundamental  to  the  successful  utilization  of  possibly  important  future  technologies  such  as 
frequency  domain  optical  storage.  7 


Inhomogencously  broadened  absorption  lines  arc  usually  regarded  as  smooth,  Gaussian 
£ 

profiles.  Recently,  we  demonstrated  the  interesting  fact  that  significant  fine  structure  is  a 
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fundamental  property  of  such  lines,  i.e.  inhomogeneous  lines  are  not  smooth  Gaussian 


functions  if  examined  closely.  The  variations  from  the  Gaussian  shape  of  the  overall 


absorption  coefficient  observed  in  a  given  frequency  interval  Av  are  due  to  statistical  changes 
(number  fluctuations)  in  the  absorber  spectral  density  over  Av  .  The  term  "statistical  fine 


structure"  (SFS)  has  been  coined  to  refer  to  this  static,  time-independent, 


frequency-dependent  structure.  It  has  been  shown  that  SFS  can  be  detected  in  a  high 


concentration  sample  by  using  a  zero-background  technique,  namely  laser 


frequency-modulation  (FM)  spectroscopy.  In  this  work  we  will  describe  the  theoretical  and 


experimental  properties  of  SFS  at  a  fixed  temperature,  1.4  K.  We  show  how  the  value  of 
the  underlying  homogeneous  linewidth  and  information  about  the  probability  distribution  of 


center  frequencies  can  be  derived  from  the  SFS  spectra  using  autocorrelation  analysis,  a 


procedure  which  is  fundamentally  different  from  hole-burning  or  coherent  transient 


techniques. 


II.  Theory  of  Statistical  Fine  Structure  in  Inhomogeneously  Broadened  Lines 


A.  Central  Concepts  and  Definitions 


The  appearance  of  statistical  fine  structure  on  spectral  lines  is  intimately  connected  with 


the  microscopic  interactions  responsible  for  IB.  It  is  generally  accepted  that  localized 


perturbations  in  the  vicinity  of  each  absorbing  center  such  as  random  strains,  electric  fields, 


and  other  perturbations  due  to  defects  in  the  host  matrix  arc  responsible  for  the  appearance 


of  IB.  In  a  reasonable  approximation  these  perturbations  shift  and/or  split  the  absorption 


line  for  each  center  without  introducing  specific  broadening.  Thus  it  is  the  distribution  of 


center  frequencies  of  absorption  that  controls  the  actual  shape  of  the  inhomogeneous  line. 


F.ach  center  in  the  probed  volume  is  assumed  to  contribute  a  single  l.orcntz.ian  profile 


to  the  total  absorption  coefficient.  I  ct  I  (v)  =  (y/2zr)/[v2  t-  (y/2FJ  signify  a  single  (normalized) 
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Lorentzian  absorption  of  homogeneous  linewidth  y,  the  full  width  at  half-maximum 
(FWHM).  We  restrict  without  loss  of  the  essential  physics  to  the  case  of  IB  for  a  single  site, 
using  the  standard  definition  of  a  single  site  as  a  well-defined  orientation  of  the  absorber  in 
the  host  matrix.  If  a  given  center  in  the  site  acquires  a  center  frequency  of  absorption  v0 
due  to  the  perturbing  fields  in  its  vicinity,  its  absorption  function  is  Lv(v  -  v0).  The  final 
shape  of  the  complete  inhomogeneous  line  profile  is  due  to  the  way  in  which  the  various 
absorbers  take  on  specific  center  frequencies  for  absorption  according  to  the  microscopic 
values  of  the  local  perturbing  fields.  To  be  precise,  we  define  the  spectral  density  of  absorbers 
in  a  fixed  probe  volume  G(v)  by  requiring  that  G(v)dv  be  the  number  of  absorbers  with  center 
frequencies  in  the  small  interval  dv  at  v.  In  previous  work,  G(v)  has  generally  been  assumed 
to  be  a  smooth  Gaussian  function  with  FWHM  r  ,  usually  called  the  inhomogeneous 
linewidth.  In  many  IB  lines,  the  ratio  T/y  can  be  I03  or  more.  In  this  paper,  we  examine 
the  implications  of  our  recent  experiments  showing  that  G(i>)  in  fact  contains  significant, 
fundamental  fine  structure. 

The  total  (inhomogcncously  broadened)  absorption  coefficient  may  be  written 

r  +oo 

a(v)=s  G(x)Ly(v  -  x)dx  =  s  [G  *  I.y](v),  (1) 

-oo 

where  s  is  the  integrated  absorption  strength  per  molecule  and  the  asterisk  signifies 
convolution.  Wc  define  the  average  number  oT  centers  per  homogeneous  linewidth,  NH,  by 
NM  =  <  f  G(v)dv  >  ,  where  the  integration  is  performed  over  the  homogeneous  linewidth  y 
and  the  brackets  mean  to  compute  the  expectation  value,  i.c.,  the  average  value  of  the 
integral  over  many  equivalent  spectral  ranges  of  width  y  .  N,,  is  often  useful  for  estimating 
purposes  because  spectral  structure  on  G(v)  much  narrower  than  y  cannot  be  deduced  front 
measurements  of  afv)  due  to  the  convolution  operation  in  Fqn.  (1). 


In  this  paper,  the  various  absorbing  centers  in  the  sample  volume  are  assumed  to  be 
independent  and  isolated.  When  NJ(  ^  1,  it  is  fairly  obvious  that  fine  structure  would  be 
clearly  evident  on  a(v).  However,  this  paper  is  primarily  concerned  with  situations  in  which 
N„^>I  so  that  the  total  absorption  at  each  frequency  in  the  inhomogeneous  line  is  due  to 
the  superposition  of  many  absorbing  centers.  We  will  show  that  even  in  this  "high 
concentration,  but  still  isolated  center"  regime,  significant,  detectable  SFS  remains  and  does 
not  average  away. 

The  fine  structure  on  the  spectral  density  G(v)  may  arise  from  one  of  two  sources:  the 
presence  of  special  center  frequencies  that  are  more  probable  ("microsites")  or  from 
fundamental  statistical  variations  in  the  total  absorption  due  simply  to  random  selection  and 
superposition  (sometimes  called  number  fluctuations).  To  quantify  this,  we  define  the 
probability  distribution  function  (pdf)  for  center  frequencies,  p(v),  by  requiring  that  p(v)dv 
be  the  probability  that  a  randomly  selected  absorber  in  the  fixed  probe  volume  will  have 
center  frequency  in  the  range  dv  at  v.  This  pdf  is  a  direct  result  of  the  large  variety  of  local 
environments  and  center  frequencies  that  arc  accessible  to  a  randomly  chosen  center  in  the 
probe  volume.  We  distinguish  two  possibilities  as  shown  schematically  in  Figure  1.  In  case 
(a),  p(v)  is  a  smooth  Gaussian  profile  indicative  of  a  normally-distributed  set  of  local 
perturbations  for  the  single  site.  In  case  (b),  p(v)  contains  departures  from  the  normal 
distribution  implying  that  certain  center  frequencies  arc  more  probable  than  others  due  to 
the  presence  of  special  defects  or  "micrositcs".  This  portrayal  of  p(v)  for  microsites  is  only 
schematic  in  that  the  actual  departures  from  a  Gaussian  function  could  be  cither  broader 
or  narrower  than  shown  in  the  figure. 

A  central  goal  of  this  paper  is  to  determine  whether  our  observations  of  SI'S  give  support 
to  one  or  the  other  form  for  p(v)  shown  in  figure  1.  In  doing  so,  it  is  important  to  realize 
that  the  actual  absorber  spectral  density  G(v)  measured  for  a  fixed  probe  volume  is  just  a 
statistical  sample  of  the  underlying  pdf,  p(v).  This  means  that  if  the  absorbing  centers  could 


be  removed  from  the  host  and  re-doped  into  the  sample  many  times,  the  resulting  G(v) 
functions  measured  for  the  same  volume  from  the  various  samples  would  give  information 
about  the  form  of  p(v). 

B.  Scaling  of  Statistical  Fine  Structure 

It  is  instructive  to  simulate  the  form  of  the  fine  structure  that  would  be  expected  due  to 
statistics  alone  for  p(v)  given  by  Fig.  1(a).  Consider  a  fixed  frequency  interval  Av  within  the 
inhomogeneous  line  that  satisfies  r^>Av$>y.  On  the  interval  Av,  p(v)  would  be  constant  and 
frequency  independent.  Since  the  various  absorbing  centers  arc  independent  and  isolated,  a 
uniform  random  number  generator  can  be  used  to  select  center  frequencies  v0,  and  the  total 
absorbance  can  be  generated  by  adding  corresponding  Lorentzian  functions  Fy(v  -  v0)  . 
Figure  2,  trace  (a),  shows  a  simulated  a(v)  for  only  two  absorbing  centers  each  with 
y  =  0.IAv  corresponding  to  TM,,  =  0.2. 

Traces  2(b)-(c)  show  a  simulation  of  a(v)  for  N,,  =  1,  10,  100,  and  104,  respectively.  Let 
Aa  and  a  signify  the  rms  amplitude  and  mean  value  of  a  on  the  interval  Av.  The  figure  shows 
that  a  grows  linearly  with  NM,  while  the  relative  rms  deviations  Aa/a  fall  approximately  as 
(N||)",/2.  This  is  reminiscent  of  the  well-known  result  from  statistical  physics  10  that  the 
relative  fluctuation  of  an  additive  physical  quantity  over  statistically  independent  systems  falls 
as  the  inverse  square  root  of  the  number  of  systems. 

To  sec  this  more  precisely,  allow  the  interval  Av  to  be  divided  up  into  smaller  intervals 

of  width  y  and  compute  N,,  =  f  G(v)dv  for  each  interval.  Since  the  various  intervals  arc 

(v) 

approximately  independent  and  the  overall  pdf  is  uniform,  the  values  of  N,,  can  be  viewed 
as  successive  samples  of  an  additive  random  variable;  hence  (AN,,)/N,,  ^  (Nlt)  12  where 
ANn  is  the  rms  fluctuation  of  N„  .  Since  on  each  interval  the  absorption  coefficient  is 
directly  proportional  to  the  number  of  centers,  we  conclude  that 

Aa/a  =  (ANn)/N||  s  (NM)  1/2  .  Since  these  deviations  of  the  absorption  coefficient  from  its 
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average  value  follow  only  from  the  statistics  of  formation  of  the  line,  we  call  the  resulting 
frequency-dependent  deviations  "statistical  fine  structure". 

Methods  of  detecting  SFS  that  are  sensitive  to  the  overall  absorption  such  as  fluorescence 
and  transmission  spectroscopy  would  find  the  largest  relative  fluctuations  at  small  Nn. 
However,  detecting  a  very  small  number  of  centers  in  the  presence  of  interfering  background 
from  the  host  matrix  is  very  difficult.  "  On  the  other  hand,  if  detection  techniques  are  used 
that  measure  only  A  a,  i.e.,  the  deviations  of  a  from  its  average  value,  the  signal  would  grow 
as  (Nu)1/2,  and  easy-to-manufacture  high  concentration  samples  would  yield  more  easily 
detectable  SFS  (as  long  as  the  concentration  is  not  so  large  that  no  light  reaches  the 
detector).  To  show  that  the  absolute  signal  in  the  NH^>1  regime  is  clearly  larger,  Figure  1(0 
shows  the  absorption  coefficient  of  1(e)  with  the  background  subtracted.  The  pcak-to-peak 
variations  in  a  are  1,600  units  compared  to  only  6  units  for  trace  1(a). 

Assuming  that  a  zero-background  detection  technique  is  used  that  measures  only  (A a)l 
where  /  is  the  sample  length,  one  may  easily  compute  the  dependence  of  the  signal  on  the 
sample  dimensions  and  center  concentration: 


(Aa)/~  a(Nn)  'Z/A  =  <r(p//A)  , 


where  a  is  the  peak  absorption  cross  section  per  center,  A  is  the  beam  area,  and  the  volume 
density  per  homogeneous  linewidth  is  p  =  Nt(/A/  .  Therefore,  (Aa)/  and  hence  the  SFS  signal 
increase  if  the  concentration  of  absorbers  or  the  sample  length  increases,  and  increase  for 
smaller  laser  spots.  Further,  centers  with  higher  cross  section  lead  to  larger  SFS  signals. 

For  a  given  beam  area,  the  maximum  signal  occurs  when  the  sample  length  /  is  equal  to 
the  Rayleigh  range.  In  this  limit,  (Aa)l  —  n(2p/Xy12,  where  X  is  the  the  optical  wavelength, 
and  the  signal  is  independent  of  the  focal  diameter.  In  this  regime,  the  probing  intensity  can 
be  decreased  without  loss  of  signal  by  choosing  a  larger  beam  area  and  a  correspondingly 
longer  sample.  Of  course,  an  ideal  sample  geometry  for  detection  of  SI'S  would  be  a  doped 
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optical  fiber  in  which  a  small  focused  spot  can  be  maintained  over  an  arbitrarily  long 
distance.  For  more  conventional  samples,  the  sample  length  is  limited,  and  small  values  of 
A  are  required  to  observe  SFS. 

C.  Detection  Using  Laser  FM  Spectroscopy 

9  12 

In  this  paper,  SFS  on  a(v)  will  be  detected  using  laser  FM  spectroscopy.  ’  The  specific 
experimental  apparatus  will  be  described  in  Section  III;  here  we  only  describe  how  the  use 
of  this  technique  enables  us  to  detect  the  structure  on  a(v).  In  FM  spectroscopy,  a 
phase-modulated  light  beam  consisting  of  a  carrier  at  the  lascr  frequency  v  and  two  sidebands 
at  (v  +  vm)  and  (v  —  vm)  is  generated  by  passing  the  lascr  beam  through  an  electro-optic  phase 
modulator  driven  at  vm.  The  modulating  frequency  vm  is  chosen  to  be  larger  than  several 
MHz,  so  that  detection  may  be  performed  at  frequencies  where  the  lascr  noise  consists 
essentially  only  of  quantum  shot  noise.  13  Such  a  phase-modulated  light  beam  contains  no 

1 4 

amplitude  modulation  as  long  as  residual  AM  effects  are  not  present.  The 
phase-modulated  beam  passes  through  the  sample,  and  amplitude  modulation  then  appears 
whenever  the  amplitudes  and/or  phases  of  the  two  sidebands  arc  disturbed  by  a 
frequency-dependent  spectral  feature  in  the  sample.  A  high-speed  detector  and 
phase-sensitive  electronics  are  used  to  detect  any  amplitude  modulation  on  the  transmitted 
beam.  If  there  is  no  spectral  feature  present,  the  photocurrent  is  not  modulated;  thus  the 
FM  signals  appear  on  zero  background. 

A  crucial  feature  of  the  FM  technique  is  that  by  simple  adjustment  of  an  electronic  phase 
shifter,  detection  may  be  performed  in  cither  quadrature  with  respect  to  the  modulating 
signal.  Adjustment  of  the  detection  phase  controls  whether  the  detected  signal  is  proportional 
to  the  absorptive  or  the  dispersive  effects  of  the  underlying  spectral  feature,  flic  detected 
signal  varying  as  cos  (27rvmt)  ,  F,(v),  is  proportional  to 
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where  Pc  is  the  laser  power  on  the  sample,  M  is  the  modulation  index  measuring  the  strength 
of  the  sidebands  relative  to  the  carrier,  /  is  the  sample  length,  and  a0  is  the  background  value 
of  a.  This  expression  is  strictly  correct  in  the  limit  of  small  M;  higher  order  expressions 
exist  5  but  are  not  important  for  the  purposes  of  this  paper.  Thus  the  cosine  phase  signal 
F[(v)  measures  the  difference  in  a/  at  the  two  sideband  frequencies,  as  long  as  <xj  is  not  too 
large  (i.e.,  as  long  as  the  exponential  factor  in  Bqn.  (3)  is  near  unity). 

The  component  of  the  detected  photocurrent  varying  as  sin(2^vmt),  I72(v),  is  proportional 


F20)~MPoc  [4>(v  +  vm)  +  4>(v  -  vm)  -  20(v)],  (4) 

where  4>(v)  =  n(v)/(27tv/c),  n(v)  is  the  index  of  refraction,  and  c  is  the  speed  of  light.  Thus, 
F2  measures  the  second  difference  of  the  frequency-dependent  index  of  refraction  at  the 
frequency  of  the  carrier  and  the  two  sidebands,  as  long  as  a0l  <  1  . 

For  either  choice  of  phase,  the  detected  FM  signal  approaches  one  of  two  well-defined 
limits  depending  upon  the  size  of  the  modulating  frequency  vm  compared  to  the  width  of  the 
spectral  feature.  When  vm  >  y,  it  is  easy  to  see  that  only  one  sideband  or  the  carrier  interacts 
with  the  spectral  feature  at  a  time.  Figure  3(a)  and  (b)  shows  the  resulting  FV1  signals  for 
a  single  Lorcntzian  profile  Ty  in  both  phases.  In  the  cosine  phase  (Fig.  3(a)),  the  FM  signal 
consists  of  a  pair  of  copies  of  the  original  I.y(v),  one  positive  and  one  negative,  separated 
by  2vm.  In  the  sine  phase  (Fig.  3(b)),  three  copies  of  the  dispersion  due  to  I.v(v)  arc  generated 
as  shown,  spaced  by  vm.  Further,  when  the  modulating  frequency  is  much  larger  than  y,  the 
FM  signals  have  an  amplitude  that  is  independent  of  the  value  of  vm.  In  the  opposite  limit 
of  small  modulating  frequency,  traces  3(c)  and  3(d)  show  that  for  cither  phase  the  resulting 
I'M  signals  approach  the  derivative  of  the  Lorcntzian,  and  in  fact  the  amplitude  of  the  signals 
goes  to  zero  linearly  with  vm.  More  examples  of  1  M  lineshapes  for  intermediate  values  of 
vm/y  arc  shown  in  Ref.  (12). 


Since  the  largest  FM  signals  occur  for  vmJ>y,  this  regime  of  FM  spectroscopy  was  used 
whenever  possible  to  detect  SFS.  One  may  think  of  Fqns.  (3)  and  (4)  as  defining  an  "I'M 


transformation"  that  is  applied  to  the  total  absorption  in  Fqn.  (I). 

D.  Autocorrelation  of  SFS  Spectra 

Estimates  of  the  homogeneous  linewidth  y  can  be  extracted  from  the  FM -detected  SFS 
spectra  by  exploiting  the  properties  of  autocorrelation  functions.  This  occurs  because  the 
total  SFS  signal  is  a  linear  superposition  of  many  single  Lorcntzian  FM  spectra  with  different 
center  frequencies.  In  a  sense,  the  autocorrelation  operation  removes  the  information  about 
the  individual  center  frequencies  in  order  to  reveal  the  underlying  homogeneous  width. 
Consider  the  autocorrelation  of  Eqn.  (1): 

P  +oo 

[a*a](v)  =  a(v')a(v'  +  v)dv',  (5) 

— OO 

where  the  star  signifies  correlation  as  distinct  from  convolution.  A  useful  identity'6  for 
arbitrary  functions  f,  g,  h,  and  j  is 

(f*g)*(h  *j)  =  (f*h)  *(g*j).  (6) 

Further,  the  autocorrelation  of  a  Lorcntzian  with  itself  is  another  Lorcntzian  of  twice  the 
width: 
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[Ly*Ly](v)=  I.2>.(v). 

The  result  of  applying  these  relations  to  Fqn.  (5)  is 

O*or](r)  =  s2  [(Ci*G)  *  L2y](r) 


(7) 
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Since  the  absorber  spectral  density  G  can  contain  frequency-dependent  structure,  the 
autocorrelation  G*G  is  not  simple  in  general.  Progress  can  be  made,  however,  by 
considering  the  expectation  value  of  Fqn.  (8), 


<  aka.  >  (v)  =  s2  [  <  G*G  >  *  L2  ](v), 


where  the  expectation  value  is  taken  over  many  spectra  a(v)  derived  from  the  same  pdf. 
Under  the  crucial  assumption  that  the  pdf  is  uniform  (Figure  1(a)),  <  G*G  >  approaches 
the  autocorrelation  of  a  broad  Gaussian  function.  The  result  near  v  =  0  is 


2 

<  a*a  >  (v)  ~  s  L2y(v)  +  (  terms  much  wider  than  y), 


which  near  the  origin  is  a  single  Lorentzian  profile  with  FW1IM  equal  to  2 y  . 

liquation  (10)  applies  for  the  actual  absorption,  whereas  the  measured  SFS  spectra  are 
related  to  a(v)  by  the  FM  transformation,  Fqns.  (3)  or  (4).  It  would  be  ideal  to  apply  an 
"inverse  FM  transform"  to  the  FM  data  in  order  to  study  the  actual  shape  of  a(v).  Such  a 
transformation  can  be  derived  by  Fourier  techniques  for  the  cosine  phase  with  the  result 

a(v)  ~  $r-,[F,(y)/2i  sin(2jrvmy)],  (II) 

where  9>~x  stands  for  the  inverse  Fourier  transform  and  F,(y)  =  ^[F,(v)].  Unfortunately, 
Fqn.  (11)  is  difficult  to  compute  for  actual  SFS  spectra  acquired  over  a  limited  spectral  range 
Av  due  to  the  lack  of  knowledge  of  a(v)  on  the  boundaries  of  Av  and  the  presence  of  multiple 
poles  in  the  inverse  Fourier  transform. 

Since  it  is  a  nontrivial  problem  to  invert  the  FM  operation  for  SFS  spectra,  we  next 
consider  the  autocorrelation  of  the  FM  cosine  phase  signal,  [!■',★  I- i](v)  .  We  further  restrict 
to  the  large  modulating  frequency  regime  vm^y  ,  although  similar  results  can  also  be  derived 
for  the  opposite  limit.  For  c  ~  1,  the  SFS  FM  signal  can  be  written 
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where  in  this  equation  and  w'hcrc  necessary  below  we  depart  from  rigor  in  the  notation  in 
order  to  explicitly  display  the  center  frequencies  of  the  I.orcntzian  functions.  Using  the  fact 
that 

r+oo 

Mv  “  vm)*I-v(v  +  vrn)  =  !-,(»'  -  vm)L,K  +  Vm  +  v)dv*  =  Lj,(v  +  2vm),  (13) 

J -or, 

the  result  is 

[^★F,](v)  =  s2  (G*G)  *  [  -  L2).(v  +  2vm)  +  2l.2y(v)  -  l.2y(v  -  2vm)].  (14) 

The  right-hand  side  of  the  convolution  operation  will  be  recognized  as  the  autocorrelation 
of  F,  for  a  single,  isolated  Forcntzian  absorption.  The  bracketed  expression  is  schematically 
shown  in  Figure  4(a)  for  positive  frequencies  only.  In  this  and  in  other  figures,  we  usually 
plot  the  normalized  form  of  Hqn.  (5)  for  convenience. 

As  was  the  case  for  a*a,  it  is  useful  to  take  the  expectation  value  of  the  autocorrclatcd 
SFS  FM  signal  in  order  to  simplify  the  influence  of  G*G  .  Under  the  assumption  of  no 
structure  on  p(v)  over  spectral  ranges  narrower  than  vm,  the  function  <  G*G  >  is  relatively 
constant  over  the  range  of  the  convolution  where  the  bracketed  term  is  large,  leaving 

<  F,*F,  >  (r)  ~  [  -  I.2y(v  +  2 vj  +  2I.2y(v)  -  I.2y(v  -  2vm)]  (15) 

as  the  dominant  form  near  the  origin  v  =  ().  Thus  an  estimate  for  y  can  be  derived  from  the 
I  W1IM  or,  more  reliably,  from  the  second  derivative  of  <  I  ,*1,  ^  at  the  origin.  In  words, 
the  expectation  saluc  ol  the  autocorrelation  of  the  F(  signal  lias  a  FWIIM  equal  to  twice 


that  for  the  underlying  homogeneous  absorption  lines.  The  approximations  made  are  most 
correct  near  v  =  0,  so  a  reasonable  formula  for  the  homogeneous  width  is 


V 


<  F,*F,  >  (0) 
<F,*F,>''(0) 


(16) 


where  the  double  prime  signifies  second  derivative. 


For  the  dispersion  signal  present  in  the  sine  phase,  we  define  the  normalized  index  change 
for  a  single  isolated  Lorentzian  at  the  origin  by 


—  vln 

y;  L,  • 


V  +  (y/2) 


(17) 


This  function  also  has  the  useful  property  that 


[Dy*Dy](v)  =  LJy(v). 


(18) 


For  the  complete  inhomogeneous  line  the  phase  shift  is  then 


<t>(v)  =  skn/  [G  *  Dy](v),  (19) 

where  k0  is  the  wavcvcctor  at  line  center.  Using  the  expression  for  the  FM  signal  in  the  sine 
phase  (Fqn.  4)  and  an  analysis  similar  to  that  for  F,  ,  we  find 

<  F2*F2  >  (v)  ~  D-2y(v  +  2vrn)  -  4F2 y(v  +  vm)  +  6L2y(v)  -  4L2y(v  -  vj  +  L2y(v  -  2vm)].  (20) 


The  positive  frequency  portion  of  the  bracketed  expression  is  shown  in  Figure  4(b),  which 
is  simply  the  autocorrelation  function  of  the  F2  spectrum  for  a  single  Lorentzian  as  shown 
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In  some  experiments,  it  is  necessary  to  consider  the  case  where  the  modulating  frequency 
is  much  less  than  the  homogeneous  width.  In  this  case  for  either  detection  phase,  the  FM 
signal  is  composed  of  a  superposition  of  the  first  derivative  of  many  Lorentzian  profiles  with 
randomly  located  center  frequencies.  The  expected  autocorrelation  of  the  FM  signals  may 
be  computed  by  taking  the  vm  <  y  limit  of  Eqn.  (14): 


<  F*F  >  (v)  ~ 


(v2 


where  F  stands  for  either  F,  or  F2. 


(21) 


Since  the  size  of  the  FM  signals  depends  upon  the  ratio  y/vm,  it  may  also  be  possible  to 
determine  the  homogeneous  width  by  measuring  the  dependence  of  the  rms  value  of  the  SFS 
spectra  on  modulating  frequency.  By  inspection  of  the  definition  of  the  autocorrelation 
function,  the  mean  square  value  of  a  FM  spectrum  is  equal  to  the  value  of  the 
(unnormalized)  autocorrelation  function  at  v  =*  0.  Using  Fqn.  (15),  then 

rms[rl(v)]'2(^-)',i—  .  ,22) 

The  value  of  y  could  then  be  determined  by  fitting  experimental  data  to  this  functional  form. 
This  method  for  determining  y  will  be  examined  in  Section  IV  along  with  the  method  arising 
from  Fqn.  (16). 

III.  Experimental 


The  low  temperature  electronic  absorption  spectrum  of  pcntaccnc  in  single  crystals  of 

17  IX 

p-tcrphcnyl  has  been  studied  in  detail.  '  At  low  concentrations,  pcntaccnc  exhibits  four 
S  «-S  site  origins  of  virtually  identical  intensity,  width  and  shape,  corresponding  to  four 
incquivalcnt  orientations  of  pcntaccnc  in  the  low  temperature  />-tcrphcnvl  lattice.  These 


pentacene  origin  lines  will  be  referred  to  as  0|  -()4 ,  with  increasing  index  denoting  increasing 
center  frequency  value.  The  absorption  line  for  each  site  has  an  overall  Gaussian  profile 
indicative  of  inhomogeneous  broadening,  with  a  T  ^  40  cm'.  In  previous  photon-echo 
studies,  the  homogeneous  linewidth  y  at  1.4  K  for  all  sites  has  been  measured  and  been 
shown  to  be  determined  solely  by  the  radiative  lifetime,  with  values  of  7.8,  7.3,  18.1  and  18.5 
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MHz  for  O,  -04 ,  respectively. 

Samples  of  pentacene  in  p-terphenyl  were  prepared  from  twice  sublimed  pentacene  and 

extensively  zone-refined  /i-terphenyl.  Dilute  single  crystals  were  grown  using  the  Bridgman 

method.  Polycrystalline  samples  were  made  by  placing  a  small  amount  of  powdered 

pentacene//?-terphenyl  mixture  between  glass  cover  slips  and  heating  on  a  hotplate  at  ~ 

200°  C  until  molten,  then  allowing  them  to  cool  fairly  rapidly.  Pentacene  concentrations 

ranged  from  1  x  10  5  to  2  x  10~7  mole/molc,  which  are  within  the  regime  where  the  molecules 

20 

may  be  regarded  as  isolated  absorbers.  Sample  thicknesses  ranged  from  10-20  ^m  for 
polycrystallinc  samples  and  250-400  nm  for  single  crystals,  yielding  low  temperature  peak 
optical  densities  between  0.01  and  0.15  for  the  pentacene  O,  line.  17 

SFS  spectra  were  acquired  using  the  experimental  apparatus  shown  in  Figure  5.  Fight 
from  a  single-frequency,  actively  stabilized,  tunable  CW  dye  laser  DF  (Coherent  599-21,  jitter 
linewidth  ^  2.8  MHz)  was  first  passed  through  an  amplitude  stabilizer  AS  (Cambridge 

Research  FS-100)  to  reduce  laser  intensity  fluctuations,  and  then  frequency-modulated  by  an 

* 

electro-optic  phase  modulator  FO  (AD  P,  Fascrmetrics  1239C).  The  modulated  beam  was 
reflected  off  the  input  mirror  of  a  confocal  Fabry-Pcrot  ctalon  FP  (Tropcl  240,  cither  1.5 
GHz  or  7.5  GHz  FSR).  When  the  ctalon  was  tuned  so  that  the  laser  scanned  over  the  dip 
in  reflectance  due  to  the  ctalon  fringe,  a  large  I'M  signal  was  produced  which  was  used  as 
a  standard  for  measuring  detection  sensitivity  and  for  adjustment  of  the  rf  phase.  For 
acquisition  of  SFS  spectra,  the  ctalon  spacing  was  adjusted  so  that  no  resonance  fringe 
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occurred  in  the  laser  scan  range  and  for  this  reason,  scans  were  limited  to  a  range  less  than 
the  FSR  of  the  ctalon.  The  laser  was  typically  set  to  scan  at  its  maximum  rate,  0.25  s/scan. 

The  laser  intensity  was  selected  using  a  combination  of  fixed  and  variable  neutral  density 

attenuators  VA,  and  then  focused  on  the  sample  S  using  a  100  mm  focal  length  achromatic 

doublet  lens  to  a  spot  ~  20  //m  in  diameter.  The  total  laser  power  on  the  sample  was 

approximately  5  n'W,  adjusted  to  compensate  for  sample  absorption,  scattering,  and  cryostat 

losses  in  order  to  maintain  1.5  /iW  at  the  detector,  l  hc  transmitted  light  was  collected  and 

focused  onto  the  detector  PD  (a  Si  avalanche  photodiode,  cither  Analog  Modules 

LNOR-275MHz-l5kV/A  or  an  RCA  C 30950 1  with  integral  preamplifier).  The  detector 

output  was  amplified  by  A,  (Q-Bit  538,  34  dB)  and  sent  to  the  rf  port  of  a  double-balanced 

rf  mixer  M  (Mini-Circuits  ZFM-3H).  The  rf  modulation  frequency  was  provided  by  a 

precision  rf  signal  generator  O  (Marconi  Instruments  20I9A)  whose  output  was  split  by  a 

two-way,  0°  power  splitter  PS  (Mini-Circuits  Zl  SC-2-5).  One  splitter  output  port  was  used 

to  drive  the  local  oscillator  port  of  the  mixer;  the  other  output  port  was  amplified  by  A2  (HNI 

503L,  40  dB)  in  scries  with  A,  (Hewlett-Packard  230B)  and  used  to  drive  the  electro-optic 

modulator,  typically  with  30  V  rms  across  a  50  £2  load.  T  his  provided  a  modulation  index 

M  ~  0.3.  Adjustment  of  the  rf  phase  of  the  detector  output  at  the  mixer  was  accomplished 

by  slightly  detuning  the  center  frequency  of  A,.  In  this  way,  the  I'M  transform  of  cither  the 

absorbance  (F,  )  or  dispersion  (F2)  of  an  SI'S  spectrum  could  be  obtained,  flic  mixer  IF 
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output  was  amplified  by  PMPA  (a  low-noise  post-mixer  pre-amplifier,  dc-20  kHz)  and 
band-pass  filtered  by  BPF  (  Tektronix  7A22  differential  amplifier  in  a  7904  oscilloscope, 
bandpass  -3dB  points  1  Hz  -  300  Hz,  typically)  before  being  digitized  and  averaged  by  DS 
(cither  a  Tektronix  468  oscilloscope,  a  TcCroy  3500SA  transient  digitizer  or  a  Data  Precision 
Data  6000  waveform  analyzer). 


It  should  be  kept  in  mind  throughout  this  work  that,  unless  specifically  noted,  no 
significance  should  be  placed  on  the  absolute  intensities  of  the  SI'S  signals  since  the  detector 
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and  overall  amplification  varied  from  one  experiment  to  another  in  order  to  maximize  the 
signal-to-noise.  This  does  not  in  any  way  affect  the  conclusions  which  are  made  herein. 

IV.  Results  and  Discussion 

As  previously  stated,  SFS  is  a  static,  time-independent  property  of  inhomogeneous 
absorption  lines  at  low  temperatures.  This  is  illustrated  for  the  case  of  the  O,  line  of 
pentacene  in  /7-terphenyl  in  Figure  6.  The  laser  center  frequency  was  repetitively  scanned 
over  a  2  GHz  frequency  range,  with  SFS  data  acquired  during  the  both  the  forward  and 
backward  scan.  In  this  figure,  an  FM  signal  of  1  V  corresponds  to  a  value  of  (Aa)/  of 
approximately  10~3.  The  large  peaks  at  0  and  2  MIIz  arc  end-of-scan  artifacts  and  will  be 
discussed  below.  The  smaller,  narrower  features  between  the  large  peaks  are  the  actual  SFS 
on  the  (),  line.  The  mirror  symmetry  between  the  up-scan  and  down-scan  regions  of  the 
spectrum  is  convincing  evidence  that  the  SI'S  peaks  are  reproducible  features  in  the  frequency 
domain;  their  static  nature  is  further  demonstrated  by  the  similarity  of  successive  averages 
taken  several  minutes  apart,  one  of  which  is  shown  in  Figure  6  offset  vertically  for  clarity. 

Figure  6  was  obtained  using  a  higher  laser  power  and  a  smaller  M  than  the  values  stated 
in  Section  III;  the  high  laser  power  produced  the  artifacts  at  0  and  2  GHz  (representing  the 
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limits  of  the  scan)  which  are  signatures  of  a  shallow  trench  caused  by  weak  hole-burning 
during  the  scanning  of  the  spectrum.  This  hole-burning  is  very  inefficient  and  the  saturation 
hole  depth  is  very  small.  For  the  results  presented  in  the  rest  of  this  paper,  the  effect  of 
hole-burning  on  the  SFS  spectra  was  made  negligible  by  reducing  the  total  laser  power  to 
1.5  mW  at  the  detector  and  increasing  the  FM  modulation  index  to  the  maximum  value  we 
could  achieve  with  our  modulator,  M  ~0.3. 

Further  information  about  the  nature  of  SFS  spectra  is  depicted  in  Figure  7,  where  the 
bottom  panel  shows  a  5  GII/  wide  scan  of  the  SI  S  near  the  peak  of  the  (),  line.  Broad  range 
SFS  spectra  like  this  arc  easily  mistaken  for  a  snapshot  of  some  noisy  signal.  However,  their 
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static,  reproducible  nature  is  demonstrated  by  the  two  traces  in  the  top  panel  of  Figure  7, 
the  lower  of  which  is  an  expansion  of  the  trace  in  the  bottom  panel;  the  upper  trace  is  a 
subsequent  acquisition  of  the  same  spectral  region  more  than  20  min  later.  The  agreement 
is  quite  good,  indicating  that  at  this  temperature  (1.4  K)  there  is  little  or  no  rearrangement 
or  annealing  of  the  perturbations  producing  the  IB  which  were  "frozen"  when  the  crystal 
was  cooled.  Most,  if  not  all,  of  the  slight  differences  between  the  spectra  are  due  to  the 
quantum  and  avalanche  noise  of  the  detection  system.  However,  some  of  the  changes  in  the 
SFS  structure  may  be  induced  by  the  probing  process  since  weak  persistent  spectral 
hole-burning  does  occur.  The  signal-to-noisc  ratio  (SNR)  of  these  and  other  averaged  spectra 
of  the  O,  SFS  in  this  work  is  ~  8-10,  measured  as  the  ratio  of  the  rms  signal  levels  when 
the  laser  is  tuned  on  and  off  the  pcntaccnc  absorption  line. 


The  possibility  of  spontaneous  rearrangement  of  the  SFS  structure  at  low  temperature 
has  not  been  ruled  out,  however.  Preliminary  experiments  to  determine  the  magnitude  of 
spontaneous  rearrangement  at  slightly  higher  temperatures  were  inconclusive  due  to  technical 
problems  associated  with  the  focused  geometry  needed  to  observe  large  SI'S  signals  in  thin 
samples.  The  slight  motion  of  the  sample  during  temperature  cycling  prevented  reliable 
measurement  of  the  same  sample  volume  before  and  after  the  temperature  excursion.  We 
did  note  however  that  the  SI  S  spectra  obtained  at  slightly  higher  temperatures  (ca.  4-10  K) 
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were  significantly  broadened.  This  is  consistent  with  the  published  dependence  of  y  on 
temperature  for  this  system.  Although  it  is  not  the  purpose  of  this  work  to  investigate  the 
temperature  dependence  of  y,  SFS  can  and  should  be  used  to  do  so  in  future  studies. 


Laser  I'M  spectroscopy  can  be  used  to  probe  cither  the  absorption  or  dispersion  of  a 
sample,  depending  on  the  relative  phases  of  the  signals  at  the  local  oscillator  (I.)  and  rf'(R) 
inputs  to  the  mixer.  When  these  signals  arc  in  phase  (the  F,  signal,  Fqn.  (3))  the  absorption 
of  the  sample  is  probed;  when  they  arc  90°  out  of  phase,  a  signal  dependent  on  the  dispersion 
of  the  sample  is  generated  (l\  ,  Fqn.  (4)).  Representative  spectra  taken  in  the  different  phases 


over  same  spectral  region  are  shown  in  Figure  8  for  the  O,  line  of  pentacene.  Traces  8(a) 
and  (b)  show  the  F,  and  F2  signals  when  vm  =  150  MHz,  respectively,  and  traces  (c)  and  (d) 
are  the  F,  and  F2  signals  when  vm  =  20MHz  .  Each  spectrum  was  obtained  twice  to  check 
the  reproducibility  of  the  SFS.  Note  that  the  "density"  of  features  in  the  F2  spectrum  in 
8(b)  is  slightly  greater  than  that  in  8(a);  this  occurs  because  the  F2  transform  produces  three 
features  from  each  Lorentzian  dispersion  feature  when  vm/y  is  large,  while  the  F,  transform 
produces  only  two  features  from  each  Lorentzian  absorption  line.  Otherwise,  the  F2  spectra 
seem  quite  similar  to  those  taken  in  F,  phase.  As  illustrated  in  8(c)  and  (d),  the  magnitude 
of  the  SFS  decreases  when  vm  decreases  because  the  FM  sideband  splitting  is  comparable  to 
the  width  of  the  features  being  probed  (see  Figure  2(c)  and  2(d)  for  comparison).  The 
detailed  dependence  of  the  SFS  signals  on  the  value  of  vm  will  be  discussed  below. 

Since  the  SFS  technique  probes  only  a  small  volume  of  the  sample,  it  is  interesting  to 
investigate  how  (or  if)  the  SFS  structure  changes  as  one  moves  transversely  to  different 
spatial  locations.  The  result  of  such  a  study  is  presented  in  Figure  9,  which  shows  a  portion 
of  the  F2  SFS  "landscape"  of  a  fixed  frequency  region  of  the  O,  line  for  a  small  range  of 
laser  spot  positions  in  a  single  crystal  of  pentacene  in  p-tcrphenyl.  The  acquisition  of  these 
data  was  accomplished  by  obtaining  101  F2  SFS  spectra,  each  in  a  20  /im  diameter  laser  spot, 
and  moving  the  position  of  the  laser  beam  on  the  sample  by  2  pm  along  a  single  spatial  axis 
between  spectra.  The  same  data  arc  presented  in  a  different  way  in  figure  10,  where  the 
SFS  intensities  arc  represented  as  various  shades  of  grev  in  a  two-dimensional  plot  of  laser 
frequency  versus  laser  spot  position,  ('lose  examination  of  the  data  reveals  that  the  features 
extend  over  several  spot  positions,  which  is  merely  a  reflection  of  the  fact  that  adjacent 
spectra  were  obtained  from  highly  overlapping  probe  volumes.  The  important  point  here  is 
that  there  arc  no  features  which  extend  over  a  long  range  along  cither  the  frequency  or 
position  dimensions.  This  is  indicative  of  the  absence  of  micrositc  structure  in  the  frequency 
range  detectable  by  the  FM  modulation  frequency  used  (58.1  Mil/).  In  other  words,  since 
the  SI'S  structure  changes  dramatically  and  apparently  randomly  from  one  sample  volume 
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to  another,  special  frequencies  that  are  more  probable  than  others  do  not  appear  to  be 
present. 

It  is  possible  to  search  for  underlying  correlations  in  the  data  of  Figs.  (9)  and  (10)  by 
computing  autocorrelation  functions  along  both  the  frequency  and  space  axes.  The  top  panel 
of  Figure  11  shows  the  average  of  101  autocorrelations  along  the  frequency  axis,  and  there 
is  no  significant  feature  in  this  autocorrelation  which  is  not  described  by  Fqn.  (20).  The 
presence  of  micrositc  structure  would  be  expected  to  produce  some  new  feature,  not  described 
by  the  theory  presented  in  Section  II.  The  bottom  panel  shows  the  average  of  512 
autocorrelations  along  the  spatial  dimension,  one  for  each  frequency  in  the  laser  scan.  This 
operation  tests  whether  a  large  peak  at  a  given  frequency  has  a  tendency  to  reappear  at  other 
spatial  locations  some  distance  away.  The  shape  of  the  lower  trace  approximates  a  (sin  x)/x 
function,  which  is  a  result  of  the  fact  that  the  SFS  structure  along  the  position  axis  has  been 
smoothed  (low-pass  filtered)  by  the  laser  spot  that  extends  over  several  spatial  positions. 

The  absence  of  sharp  peaks  in  the  lower  trace  suggests  that  spatial  locations  separated  by 
more  than  the  laser  spot  diameter  of  20  nm  are  uncorrclatcd.  Both  these  observations  give 
support  to  the  statistical  origin  for  the  observed  SFS,  i.e.,  there  is  no  evidence  to  suggest  that 
micrositcs  exist  in  this  material  over  frequency  ranges  smaller  than  approximately  100  MHz. 

We  have  also  performed  a  preliminary  examination  of  the  SFS  structure  of  the  02 , 03 
and  04  site  origins  of  pentaccne  in  /7-tcrphenyl.  Figure  12  shows  three  in-focus  1.4  K 
spectra:  trace  (a)  is  the  FM  signal  acquired  at  a  laser  frequency  far  removed  from  the 
pentaccne  absorption.  1  his  noise  level  represents  the  quantum  noise  from  the  probe  beam 
and  the  detector  avalanche  noise  at  150  MHz..  In  particular,  the  residual  electronic  noise 
with  no  light  on  the  detector  is  well  below  this  level,  and  any  residual  AM  from  the 
modulator  is  also  below  this  level  due  to  the  low  light  powers  in  use.  Trace  (b)  is  the  I'M 
SFS  signal  near  the  peak  of  the  (),  line;  and  trace  (c)  is  the  FM  SFS  spectrum  near  the  peak 
of  the  Oj  line.  The  SFS  in  the  (),  line  is  strikingly  different,  with  lower  intensity  and  broader 


spectral  features  than  for  the  O,  site,  Results  similar  to  those  obtained  from  03  occur  for 

the  04  site,  while  the  SFS  of  the  02  site  mimics  that  of  the  O,  line.  Altnough  the  measured 

homogeneous  Iinewidths  for  03  and  04  are  more  than  twice  as  large  as  those  for  O,  and 
19 

02 ,  and  therefore  should  produce  smaller  FM  signals  for  a  given  modulation  frequency, 
this  effect  alone  does  not  explain  the  far  smaller  amplitude  of  SFS  observed  on  these  sites. 
The  weaker  and  broader  SFS  for  the  two  higher  sites  may  also  be  due  in  part  to  the 
significant  increase  in  the  S,  -*  T,  intersystem  crossing  (ISC)  quantum  yields  for  the  03  and 
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04  sites  (d>ISC  ~  0.6)  ,  compared  to  that  for  the  O,  and  02  sites  (<t>ISC  a;  4x  10  3).  This 
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increase  in  d>ISC,  together  with  a  triplet  lifetime  of  70  /us  '  allows  easy  transient 
bleaching  of  the  S,  <-  S0  transition  during  probing,  thereby  significantly  reducing  the  number 
of  centers  contributing  to  the  SFS  FM  signal  and  reducing  the  amplitude  of  the  SFS 
spectrum. 

In  Section  IIB  we  argued  that  the  rms  intensity  of  the  SFS  should  be  proportional  to 
(Nu),/2  (Fqn.  (2)).  Figure  13  shows  the  rms  amplitude  of  a  large  number  of  O,  FM  SFS 
spectra  versus  N»  for  a  sequence  of  seven  polycrystallinc  thin  film  samples.  The  samples 
were  all  prepared  from  a  single  parent  concentration  of  pcntacenc  in  p-tcrphcnyl,  using 
different  sample  thicknesses  to  adjust  Nu.  The  actual  value  of  NM  for  each  sample  was 
determined  from  the  room-temperature  optical  density,  the  extinction  coefficient 
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(e  =  9.9  x  103)  and  the  low-temperature  inhomogeneous  broadening  factor  for  these 
polycrystalline  samples,  T/y  =  1.1  x  104  .  I  hc  straight  line  is  a  least-squares  fit  to  the  data 
showing  that  the  rms  amplitude  of  F,(v)  grows  as  (N,,)0 34,005  which  is  consistent  with  the 
expected  (NM)I/2  dependence.  The  vertical  scatter  in  the  points  is  understandable,  because 
each  FM  spectrum  was  acquired  over  a  different  spectral  range. 

The  rms  value  of  a  given  SFS  spectrum  also  depends  upon  the  relative  si/cs  of  vm  and 
y  (sec  Figures  3  and  8)  and,  in  principle,  it  is  possible  to  determine  the  underlying 
homogeneous  linewidth  of  an  SFS  spectrum  by  measuring  the  rms  value  as  a  function  of 


vm.  Figure  14  shows  the  measured  vm  dependence  of  the  in-focus  rms  SFS  signal  for  the  O, 
line  using  F,  .  Each  data  point  (circle)  at  a  particular  vm  represents  the  mean  rms  value  of 
a  set  of  11  in-focus  O,  SFS  spectra  taken  at  different  center  frequencies  to  provide  a  good 
statistical  sample  of  the  SFS.  This  method  reduces  any  scatter  from  the  true  rms  SI'S  value 
which  might  occur  for  a  single  spectrum  due  to  local  fluctuations  of  G(v).  The  error  bars 
represent  two  times  the  standard  deviation  of  the  set  of  individual  rms  values  obtained  at  a 
given  vm.  The  same  set  of  center  frequencies  was  used  at  each  new  vm.  The  data  are 
corrected  for  modulation-frcquency-depcndent  changes  in  the  detection  sensitivity,  most 
notably  for  changes  in  the  modulation  index.  The  two  solid  lines  in  Figure  14  arc  simulated 
curves  using  Eqn.  (22):  the  lower  line  assumes  a  value  of  y  equal  to  40  MHz,  and  the  upper 
line  is  the  analogous  simulation  using  the  previously  measured  value  of  y,  7.8  MHz.  The 
non-zero  value  of  the  simulations  at  vm  =  0  represents  the  residual  off-line  averaged 
quantum  +  avalanche  noise  level  of  13  mV.  Clearly,  the  determination  of  the  homogeneous 
linewidth  by  this  fitting  procedure  misses  the  actual  value  by  a  large  factor. 


The  departure  of  the  experimental  data  from  the  expected  7.8  MHz.  curve  can  be 
explained  in  the  following  way.  The  FM  signal  from  a  sample  is  generated  when  cither  of 
the  low  intensity  FM  sidebands  is  in  resonance  with  some  spectral  feature.  If  the  intensity 
of  the  laser  carrier  on  the  sample  is  above  the  saturation  threshold,  power  broadening  of  the 
homogeneous  linewidth  and  a  correspondingly  reduced  absorbance  in  the  immediate  vicinity 
of  the  carrier  frequency  will  occur.  If  vm  is  greater  than  the  width  of  the  transient  spectral 
hole  formed  at  the  carrier  frequency,  then  both  sidebands  will  be  probing  an  undistorted 
region  of  the  spectral  line,  as  long  as  the  transient  spectral  hole  from  the  carrier  decays  before 
cither  sideband  can  reach  the  spectral  region  of  reduced  absorption.  If  vm  is  comparable  to 
or  less  than  the  width  of  the  transient  spectral  hole  formed  at  the  carrier  frequency,  then 
both  sidebands  will  probe  the  broadened  features  near  the  carrier.  Hence  the  modulation 
frequency  dependence  at  low  vm  will  be  dominated  by  power  broadening  due  to  the  carrier, 
from  the  data  presented  in  Figure  14,  we  conclude  that  a  power-broadened  transient  hole 


tens  of  MHz  in  width  has  been  formed  around  the  carrier  frequency.  This  effect  can  be 
minimized  by  reducing  the  total  laser  intensity  so  that  the  carrier  would  no  longer  saturate 
the  absorption,  but  this  also  reduces  the  SFS  signal  amplitude  proportionally  (Eqn.  (3)).  In 
principle,  this  could  be  compensated  for  by  increasing  the  modulation  index,  but  in  these 
experiments  the  electro-optic  modulator  was  always  driven  at  the  maximum  rf  power 
available  (30  V  rms  into  50  £2).  When  the  laser  power  was  decreased,  the  SNR  was 
insufficient  to  produce  good  reproducibility. 


A  much  better  method  of  determining  y  which  does  not  necessitate  the  use  of  low 
modulation  frequencies  takes  advantage  of  the  properties  of  autocorrelation  functions  (see 
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Eqns.  (15)  and  (16)),  and  has  been  successfully  used  on  the  pcntaccnc/p-terphenyl  system. 
These  equations  relate  to  the  expectation  value  of  E,*!',;  the  autocorrelation  of  a  single  F, 
SFS  spectrum  docs  not  provide  a  large  enough  statistical  sample  of  the  SFS  structure  to  yield 
a  reliable  approximation  for  <  F,*!-',  >  and  therefore  a  reliable  value  for  y.  To  illustrate  this, 
normalized  autocorrelations  of  single  in-focus,  (>,  SFS  spectra  taken  with  vm  =  I50MIIz  for 
F,  and  F2  phases  arc  shown  in  Figure  15(a)  and  (b),  respectively.  Both  traces  show  a  large 
feature  at  the  origin  with  approximately  the  correct  FWI1M  (2y)  along  with  many  random 
fluctuations  about  the  zero  level  which  bear  little  resemblance  to  the  expected  forms  of  Eqns. 
(15)  and  (20).  However,  if  many  SFS  spectra  arc  obtained  over  different  frequency  intervals 
and  their  autocorrelations  averaged,  the  expected  dependence  becomes  quite  apparent.  Traces 
15(c)  and  (d)  show  the  results  in  the  two  phases  for  averages  of  66  and  11  spectra, 
respectively.  One  can  sec  in  these  averaged  autocorrelations  the  negative  feature  predicted 
at  2vm  for  F,  as  well  as  the  negative  feature  at  vm  and  the  positive  feature  at  2vm  for  !\ 
as  given  in  Eqns.  15  and  20.  The  number  of  spectra  needed  to  produce  an  average 
autocorrelation  which  approaches  the  proper  expectation  value  decreases  as  the  frequency 
range  of  the  scans  is  widened,  because  errors  in  the  calculation  of  the  autocorrelation  due 
to  the  finite  spectral  range  arc  minimized. 
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For  a  single  crystal  sample  with  concentration  3  x  10-7  moles/mole,  we  measure 

y  =  7.9  +  0.8  MHz  for  the  O,  site  at  1.4  K  by  analyzing  the  autocorrelations  of  six  SFS 

spectra  using  Fqn.  (16).  The  second  derivative  of  the  autocorrelation  function  at  the  origin 

2*> 

was  determined  by  using  a  standard  five-point  central  difference  formula.  Fach  spectrum 

was  acquired  over  a  different  4.8  GHz  spectral  range  using  vm  =  150  MHz.  This  experimental 
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value  for  y  is  consistent  with  the  previously  reported  value  of  7.8  +  0.6  MHz. 

Depending  upon  the  experimental  conditions,  certain  artifacts  can  distort  the  shape  of 
the  autocorrelation  function  near  the  origin  and  give  rise  to  errors  in  the  determination  of 
y  using  Fqn.  (16).  It  is  fairly  obvious  that  errors  will  occur  in  the  autocorrelation  if  too  few 
points  are  digitized  for  each  spectral  range  of  width  y.  If  the  laser  linewidth  is  not  much 
smaller  than  y,  the  absorption  profile  for  a  single  isolated  Forcntzian  will  be  the  convolution 
of  a  Forcntzian  and  the  laser  spectral  profile.  In  addition,  if  the  SNR  of  the  SFS  spectrum 
is  not  high  enough,  a  peak  at  the  origin  in  the  autocorrelation  due  to  noise  will  distort  the 
value  one  might  determine  for  the  second  derivative.  Finally,  the  spectral  range  over  which 
the  entire  SFS  spectrum  is  acquired  may  be  too  small  to  allow  approximation  of  the  integral 
in  Fqn.  (5)  by  a  finite  sum.  The  error  introduced  by  each  of  these  in  a  given  experiment 
can  be  estimated  by  extending  the  autocorrelation  analysis  in  Section  IF  In  our  experiments, 
we  determined  the  following  useful  rules  of  thumb  to  keep  the  error  in  the  use  of  Fqn.  (16) 
below  a  few  percent:  (i)  maintain  laser  linewidth  less  than  y/3,  (ii)  digitize  the  SFS  spectrum 
with  at  least  5  points  per  frequency  range  y,  (iii)  choose  vm  >  10y,  and  finally  (iv)  if  the 
inhomogeneous  line  is  large,  acquire  spectra  over  a  range  of  at  least  20vm  or  until  the  SFS 
becomes  negligible. 

V.  Conclusion 

Using  laser  I  \1  spectroscopy,  we  have  observed  SFS  for  the  origin  lines  of  the  four  sites 
of  pcntaccnc  in  />-tcrphcnyl  crystals,  and  have  presented  a  formalized  theoretical  description 


of  the  effect  and  its  analysis.  Measurement  of  the  homogeneous  linewidth  for  the  O,  site 
origin  using  autocorrelation  analysis  of  SFS  spectra  yielded  a  result  in  excellent  agreement 
with  previous  determinations.  SFS  data  taken  over  a  range  of  probe  volumes  in  a  single 
crystal  sample  to  investigate  the  underlying  probability  distribution  function  of  the 
inhomogeneous  profile  yielded  no  evidence  for  the  existence  of  microsite  structure  on  O,  over 
frequency  ranges  smaller  than  ~  100  MHz.  The  dependence  of  the  SFS  signals  on  center 
concentration,  rf  modulating  frequency  and  detection  phase  was  measured  and  found  to  be 
in  general  agreement  with  the  statistical  model  for  the  fine  structure. 

Limitations  of  the  technique  may  arise  from  the  following  facts:  (i)  the  homogeneous 

linewidth  may  be  so  large  that  high-speed  detectors  and  modulators  are  not  available  to 

achieve  the  condition  for  large  signals  vm^>y  ,  (ii)  persistent  spectral  hole-burning  may  prevent 

the  careful  detection  of  SFS  due  to  bleaching  by  the  probing  laser,  and  (iii)  the  concentration 

or  cross  section  of  the  absorber  may  be  too  small  to  yield  SFS  signals  larger  than  the  laser 

shot  noise.  Limitation  (i)  may  be  overcome  by  using  a  variation  of  I’M  spectroscopy,  the 
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two-tone  FM  technique,  which  allows  the  detection  of  much  broader  features  with  high 
sensitivity.  In  spite  of  the  remaining  technical  limitations,  SI'S  should  be  a  universal  property 
of  all  inhomogcneously  broadened  absorption  lines,  and  it  is  hoped  that  the  techniques 
described  here  will  be  applied  to  study  the  low-temperature  spectra  of  other  ions  and 
molecules  in  crystals  as  well  as  in  amorphous  media.  If  spectra  can  be  acquired  in  a  time 
short  compared  to  the  collision  time,  SFS  may  even  be  observable  on  the  Doppler-broadened 
profile  of  a  gas.  A  further  interesting  question  is  whether  or  not  SFS  can  be  detected  using 
coherent  temporal  probes  of  the  inhomogeneous  line  such  as  Tree-induction  decay. 

It  is  fairly  obvious  that  SFS  is  identical  in  origin  to  the  kind  of  spectral  structure  one 
would  expect  if  the  optical  absorption  due  to  individual  absorbers  in  a  solid  could  be 
detected.  In  more  common  samples  where  the  number  of  molecules  per  homogeneous 
linewidth  is  large,  statistical  fluctuations  of  the  spectral  density  of  absorbers  will  always  occur. 
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FIGURE  CAPTIONS 


Figure  1.  Simulated  probability  distribution  functions.  Case  (a)  is  the  smooth  Gaussian 
profile  usually  assumed  for  an  inhomogencously  broadened  set  of  absorbers;  case  (b) 
represents  a  pdf  which  includes  microsites,  or  departures  from  the  Gaussian  form  due  to 
certain  center  frequencies  having  higher  and  lower  probability  than  in  the  normal  case. 

Figure  2.  Simulated  absorption  spectra.  Each  spectrum  is  the  sum  of  single  Lorentzians  each 
with  y  =  0. 1  and  with  randomly  chosen  center  frequencies.  Traces  (a)-(e)  correspond  to  NM 
values  of  0.2,  1,  10,  100  and  104.  The  data  for  trace  (0  is  the  same  as  that  for  (e),  with  the 
scale  expanded  about  the  average  value  to  show  the  fluctuations  still  present. 

Figure  3.  Simulated  EM  spectra  from  Eq.  (3)  and  (4).  The  spectra  were  generated  from  a 
single  Lorcntzian  absorption  or  dispersion  feature  centered  at  0  MHz  with  y  =  8  Mffz. 
'Traces  (a)  and  (b)  were  obtained  with  vm  =  150  MHz  ;  traces  (c)  and  (d)  were  obtained  with 
vm=  1  MHz. 

Figure  4.  Simulated  (normalized)  autocorrelation  spectra.  These  autocorrelations  were 
calculated  from  the  simulated  EM  spectra  shown  in  Figure  3(a)  and  (b). 

Figure  5.  Experimental  apparatus.  See  text  for  description  of  symbols. 


Figure  6.  SI  S  of  the  O,  line  in  pcntaccnc  in  /Mcrphcnyl  at  1.4  K.  Each  trace  is  the  average 
of  64  scans  over  a  2  GHz  frequency  range,  flic  SI  S  is  the  smaller,  narrower  features  between 
the  large  artifact  peaks  at  0  and  2  GHz,  and  shows  mirror  symmetry  in  the  up  and  down 
scan  portions  of  the  spectrum.  The  lower  trace  is  a  second  acquisition  taken  several  minutes 
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Figure  7.  Wide  range  scan  of  SFS.  The  lower  panel  shows  the  result  of  64  averages  of  a  5 
GHz  scan  of  the  SFS  on  the  O,  line  at  1.4  K.  The  upper  panel  shows  a  portion  of  the  same 
spectrum  on  an  expanded  frequency  scale  together  with  another  acquisition  of  the  spectrum 
taken  more  than  twenty  minutes  later,  offset  for  clarity. 

Figure  8.  SFS  with  different  detection  phases  and  modulation  frequencies.  The  spectra  are 
the  average  of  64  scans,  all  over  the  same  range  of  frequencies  near  the  peak  of  the  O,  line 
at  J.4  K.  Traces  (a)  and  (h)  arc  the  F,  and  F2  responses  when  vm  =  150  MHz;  traces  (c) 
and  (d)  are  the  F,  and  F2  responses  when  vm  =  20  MHz.  Superimposed  over  each  trace  is 
another  subsequent  average  of  the  spectrum. 

Figure  9.  Surface  plot  of  SFS  spectra  at  different  laser  spot  positions.  Fach  SFS  spectrum 
is  the  average  of  64  scans  over  the  same  range  of  frequencies.  A  total  of  101  spectra  were 
obtained  in  the  Is  phase  at  vm  —  58.1  MHz  ,  each  with  the  laser  focused  to  a  20  /rm  spot 
on  the  sample  and  with  the  laser  spot  position  displaced  by  2  /an  along  a  single  spatial  axis 
between  spectra. 

Figure  10.  A  2-dimcnsional  grey  scale  plot  of  the  same  data  as  that  depicted  in  Figure  9. 
The  shades  from  white  to  black  indicate  the  value  of  the  SI'S  signal,  with  lighter  shades 
indicating  higher  values. 

Figure  11.  Normalized  average  autocorrelations  of  the  data  of  Figures  9  and  10  along  the 
two  independent  axes,  flic  top  panel  shows  the  normalized  average  autocorrelation  along 
the  laser  frequency  axis  of  all  101  SI  S  spectra.  The  bottom  panel  shows  t he  normalized 
average  autocorrelation  along  the  position  axis  at  512  discrete  laser  frequencies. 


Figure  12.  SFS  in  the  different  site  origins.  Trace  (a)  is  the  averaged  noise  level  obtained 
when  the  laser  is  tuned  off  all  absorbance  lines;  trace  (b)  is  SFS  near  the  peak  of  the  O,  line, 
and  trace  (c)  is  SFS  near  the  peak  of  ()3 .  All  spectra  were  averaged  over  64  laser  scans  in 
the  F|  phase  with  vm  =  150  MHz  and  constant  detection  sensitivity. 

Figure  13.  Dependence  of  SFS  rms  amplitude  on  NM.  Different  Nn  values  were  obtained 
by  using  different  thicknesses  of  sample  at  the  same  concentration.  The  straight  line  is  a 
least-squares  fit  to  the  data,  with  slope  =  0.54  ±  0.05  . 

Figure  14.  Dependence  of  rms  SFS  amplitude  on  vm.  Fach  circle  represents  the  average  of 
the  rms  values  of  1 1  F,  SFS  spectra  from  the  O,  line,  each  spectrum  taken  over  a  different 
frequency  range.  The  error  bars  represent  ±  two  standard  deviations  for  the  rms  values  at 
that  frequency.  The  two  solid  lines  are  simulated  curves  according  to  Fqn.  (22)  using 
y  —  7.8  and  40  MHz  for  the  upper  and  lower  curves,  respectively. 

Figure  15.  Normalized  autocorrelations  in  F,  and  I  2  •  Traces  (a)  and  (b)  show  the 
autocorrelations  of  single  F,  and  I  2  O,  SFS  spectra,  respectively.  Traces  (c)  and  (d)  arc  the 
averaged  autocorrelations  of  66  and  11  SFS  spectra,  respectively,  in  the  two  detection  phases. 
All  spectra  were  obtained  near  the  peak  of  the  (),  line  with  vm  =  150  MHz. 
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